Vibration of the Dimer on Ge(OOl) Surface Excited Coherently by STM Current 
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The vibration of the dimer on Ge(001) surface with the higher vibrational number excited co- 
herently by STM current is theoretically investigated. The coherent excitation rate of the dimer 
vibration is obtained by the Hamiltonian consisting of the terms of the electron system and the 
electron-vibration coupling terms. The transformation of the local structures reported in STM ex- 
periments is shown to be driven by the dimer vibration excited coherently by the STM current. 
The sample bias voltage of STM above which the transformation to the p(2x2) structure is able 
to be observed in the experiments, is semiquantitatively reproduced by the quasi-one-dimensional 
character of the 7r*-band. We show that the excitation rate has the term not decaying on the 
distance from the STM tip for the one-dimensional band. The contrastive diffrence of the shapes 
of the transformed region depending on the sign of the bias voltage observed by the experiments 
is explained by the difference between the quasi-one-dimensional character of the 7r*-band and the 
two-dimensional character of the 7r-band on Ge(001) surface. 

PACS numbers: 68.35.Ja, 68.37.Ef, 73.40.GK, 63.20.Kr 
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On Ge(OOl) surface, neighboring atoms form the buck- 
led dimers. The orientational arrangement of the buck- 
led dimers orders into the c(4x2) structure through the 
order-disorder phase transition about 315K Q, U0, 0, 
000 Recently, the phase manipulation of Ge(001) 
surface on Sb-doped substrate by scanning tunneling mi- 
croscopy (STM) apparatus has been reported 0, . The 
local surface structure on Ge(OOl) surface is modified re- 
versibly and hysteretically between c(4x2) and p(2x2) 
structures at 80K by the sample bias voltage V s of STM. 
The transformation of the local structure from c(4x2) 
to p(2x2) is observed at positive sample bias voltages of 
0.8V < V s and the reverse transformation frornp(2x2) 
to c(4x2) is done at negative V s < -0.7V The 
transformation takes place occasionally during the scan 
of the surface by STM. The transformation rate to the 
p(2x2) structure is reported to be proportional to the 
tunneling current of STM (Istm) 13, 13- The occasional 
transformation to the p(2x2) structure during the scan is 
restricted to the dimer row including the dimer beneath 
the tip of STM, and takes place simultaneously in the 
long region in the dimer row. In most cases, the trans- 
formed p(2x2) region in the host structure of c(4x2) con- 
tinues to the outside of the scanning area along the dimer 
row. The reverse transformation to c(4x2) presents a re- 
markable contrast to the transformation to p(2x2). The 
transformed c(4x2) region expands over the consecutive 
dimer rows beneath the tip. The region, however, is re- 
stricted within the range of a few nm away from the dimer 
beneath the tip for both directions of along and perpen- 
dicular to the dimer row. 

In the present letter, we investigate theoretically the 
vibration of the dimer on Ge(001) surface excited coher- 
ently by STM current, using the Hamiltonian consisting 
of the term of the electron system describing the tunnel- 
ing current of electrons and the electron-vibration cou- 
pling term. The excitation rate of the vibration of the 



dimer beneath the tip to the higher vibrational number 
n by the coherent process is shown to be proportional to 
Istm- We report that the excitation rate of the vibrar- 
tion of the dimer away from the tip decays much slowly 
on the distance from the tip for the direction along the 
dimer row in the case of V s > 0. The range of the dis- 
tance that the transformation is effective, is expected to 
be of the order of the mean free path of electron in the 
7r*-band. 

Previously, we have investigated the vibration of the 
dimer excited by STM current on Si(001) surface of the 
p-type substrate, where the case of negative V s has been 
treated 

G3- 

The vibrartion of the dimer beneath the 
tip on Si(001) surface is excited to the higher vibrational 
number n by the STM current through the incoherent 
ladder climbing of the vibrational states at low tempera- 
tures. The effective temperature of the vibration on the 
dimer beneath the tip strongly deviates from the sub- 
strate temperature T and reaches a few hundred Kelvin 
for the typical values of STM current at low temperatures 
lower than 20K on the boron doped (B-doped) Si(001) 
substrate. The high effective temperatures at low T cause 
the rapid flip-flop motion of the dimer on Si(001) sub- 
strate. At these low temperatures lower than 20K, the 
excitation rate of dimer vibration with the higher vibra- 
tional number n depends on the Istm in highly nonlinear 
manner on B-doped Si(001) substrate [hJ. When T in- 
creases around 20K, the deexcitation of the vibration by 
the transition from the vibrational number n to n — 1 
through the inner-band excitation of the electron-hole 
pair creation in the 7r-band becomes prominent and the 
rate of the deexcitation increase steeply in T 10]. The 
high increase of the deexcitation rate causes the steep de- 
crease of the probability of the incoherent ladder climbing 
of the vibrational states and that of the effective temper- 
ature of the dimer beneath the tip on Si(001) substrate. 
At temperatures higher than about 35 K, the incoherent 
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ladder climbing excited by STM current is not dominant 
on the excitation to the higher vibrational number n |ld| . 
We have shown that the temperature where the inner- 
band excitation becomes prominent is determined by the 
absolute value of the energy of the 7r-band top measured 
from the Fermi level of the electron system [lOj . On the 
B-doped Si(OOl) substrate, the absolute value of the en- 
ergy is about 22.5meV. 

On Sb-doped Ge(OOl) substrate used by the experi- 
ment the absolute value of the energy of the ir*- 
band bottom measured from the Fermi level of the elec- 
tron system is about 4.8meV. From this small value of 
the energy, the rate of the deexcitation of the vibration 
from n to n — 1 through the inner-band excitation in 
the 7r*-band on the Sb-doped Ge(OOl) substrate is ex- 
pected to become prominent at temperature much lower 
than that on the B-doped Si(OOl) substrate. The rate 
is estimated to be effective at temperature of the order 
of (4.8/22.5) X20K—4K. From this estimation of temper- 
ature, we know that the incoherent ladder climbing of 
the vibrational states excited by the STM current on Sb- 
doped Ge(OOl) substrate is suppressed at higer tempera- 
tures, and becomes negligible already at 8 K. Therefore, 
the transformation of the local structure observed by the 
experiment at 80K and 10K on the Sb-doped Ge(OOl) 
substrate |a, l2j can not be ascribed the flip-flop motion 
of the dimer through the incoherent ladder climbing of 
the vibrational states excited by the STM current. 

We introduce the Hamiltonian H = H c + H cv repre- 
senting the coupling between the electronic system and 
the localized vibrational system, 
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with e a = \lk\ 2 £k and T p = Tj p , where c p is the 
annihilation operators for the electronic states in the 
conduction band of the STM tip, and Ck is the anni- 
hilation operator for the electronic states in the 7r-band 
for V s < or in the 7r*-band for V s > 0. The opera- 
tors of ao = J2k 7* c fe an d Or = J2k exp(— i/c • r)~f*Ck are 
the annihilation operators for the electronic states |ao) 
spatially localized on the dimer beneath the tip and the 
electronic states \a r ) spatially localized on the dimer at 
the position r (7' 7^ 0) measured from the dimer beneath 



the tip, respectively. The operators of bo and b r repre- 
sent the annihilation operators for the vibrational state 
in the rocking mode of the dimer beneath the tip and 
of the dimer at the position r, respectively. We assume 
that the vibrational states on the dimers are described 
as the harmonic oscillator. The energy of the vibrational 
state froj is roughly estimated as 13meV 0. The diago- 
nalized term of ^ fe Skc\,Ck is rewritten by ao and c a as 
the first three terms in the third line of eq. (1), where 
c a is the annihilation operators for the electronic states 
orthogonalized to \ao). The term of Se represents the 
electron-vibration coupling. The vibration is excited or 
deexcited through the electron-vibration coupling by the 
STM current. The essentially same Hamiltonian as eq. 
(1) has been used in the previous study for Si(OOl) sur- 
face ^3 an d other surfaces 0, H2I H3I IT3 . HH Ha. HtI fli| . 

In the present study, we discuss the transition rates 
of the vibration mainly in the case of V s > 0. The 
inter-band transition of the electrons from the conduc- 
tion band of the tip to the 7r*-band can excite coherently 
the vibrational state of the dimer through the electron- 
vibration coupling. In the coherent transition, the vibra- 
tional state is excited quantumly to the excited state of 
the higher quantum number through the t-matrix. The 
rate of the coherent vibrational excitation is proportional 
to the electronic current Istm 0, • The coherent ex- 
citation rate <7o,o^n m the dimer beneath the tip from the 
ground state of the vibrational number to the excited 
state of n coupled with the inelastic electrons transition 
is obtained in the lowest order in \5e\ as 

= 2y 5^(1 - /(£«/ - £ F ))/(V - eV s - e F ) 
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x(eV s — nfiuj — E^'), 
Ge-v(e) = (e + M + - (H c + H cv ))-\ 
Goo(e) = (a o \{£ + i0+ -H^lao) 
= (e-e a - A(e) + iA(e))~ 1 , 



(2) 



where f(e) = (exp(/fe) + 1) \ G 00 {e) = (e - e a + 
iA s (e)) _1 and £f is the Fermi level of the electronic sys- 
tem in the surface. Goo an d G G - V are the Green func- 
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tions for the electronic system and for the whole sys- 
tem, respectively. A(e) is n~ 1 P J^de' A(e') / (e — e') 
where P denotes the Cauchy principal value. A(e) = 
A s (e) + A t (e) is the width of projected density of states 
for \a ). A s (e) = n £ Q \T a \ 2 5{e - e a ) and A t (e) = 
ttJ2 p |r p | 2 ^(£ — e P ) are the components of A(e) in the 
surface and in the tip, respectively. E^* is the en- 
ergy level of the bottom of the 7r*-band, and D^* is 
the width of the 7r*-band. Within the 7r*-band, A(e) is 
approximated as A(e) as A s (e), because |r Q | is to be 
much larger than |r p | in the STM experiments. \a') 
and \p') are the stationary states of H c connected to 
| a) and \p), respectively. Namely, the stationary states 
are given as \a') — \a) + (e a + i0 + — H c )~ 1 T a \a ) and 
\p') = \p) + (e P + «0+ - i/ c ) _1 r p |ao}. The broadness 
of the Fermi distribution around ep in the tt* band and 
around e-p + V s in the tip is assumed to be negligible in 
the evaluation of the integration in eq. (2). The evalu- 
ation point for the integration in eq. (2) are defined at 
e = eV s — nhuj — X n . At(e) is assumed to take constant 
value of At in the evaluations in eq. (2). Because the 
height of the energy barrier for the flip-flop motion of the 
dimer in the c(4x2) ordered structure and the height for 
the type-P defect are estimated from the model potential 
of Ge(OOl) as about 4.2 x 10 2 meV and 3.1 x 10 2 meV, 
respectively, n is estimated to be of the order of 20 ~ 30 
for the transformations. 

The surface localized states of the 7r-band and the ir*- 
band are formed mainly through the hybridization of the 
dangling bonds on the dimers. The first-principles calcu- 
lations (FPC) [i^l and scanning tunneling spectroscopy 
(STS) |1 @ show that the 7r-band and the ?r*-band 
have the two-dimensional character and the quasi-one- 
dimensional character, respectively, and the band widths 
of both bands are obtained to be about 1 eV. In quasi- 
one-dimensional bands, the density of states has large 
values near the band edges, and becomes small around 
the central region of the band. In the present study, we 
assume e a of 7r*-band to be the center of the band in 
energy: e a = E v * + D^* /2. Corresponding to the above 
mentioned character of the quasi-one-dimcnsional band, 
the energy dependence of A s (s) is simply modeled as 
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A s , (\e-£ a \ < > 
A s (e) = { A L , (0 < \e-e Q 
0, otherwise 



(3) 



We know by the STS results [ali that A L /A S and cj> are 
of the order of 2 and 50meV, respectively. Because the 
value of Al is roughly estimated about 500meV from the 
width of the 7r*-band, <j) is expected to be much smaller 
than As- 

In this model, 00,0— >n are obtained from eq. (2) for 
E^* < eV s — nhio as 



00,0- 



Tin 



n l ^ 



x (eV s - nfwj - E^* ) , (4) 

where Lj(x) — ((eVs—jTiui^ + x 2 )^ 1 . The integers t\ and 
£2 are defined by (ev s — iituv) < <f> < (ev s — (ii — l)Tiuj) and 
{eVa—lztuj < —<\> < (ev s — (£2 — l)Titjj), respectively, Ylu ^ s 
assigned to 1 for i < 0, and ev s — eV s — X n — e a . We take 
the scope to the case of ev s ~ 0, hereafter. As mentioned 
above, we assume from STS results H,0 that <C A s < 
Al- In the condition of ev s ~ 0, the voltage dependent 
terms of ev s — jTito appearing in the denominators of the 
products in eq. (4) are approximately neglected. The 
inelestic rate of ao,o-m are classified by the relative values 
of ev, to —d> or 6 as 



4n!A, 



F(ev s ) 



ihuj — E 7r *)F(ev s ), 



TrhA h 

1 (ev s < -4>) 
(A L /A s f ev ° + ^ ev s <cfy .(5) 

(A L /A S ) 4 ^ (0 < ev s ) 



3=0 



Eq. (5) shows that cr 0i o— n increases gradually with the 
powers of the bias voltage for ev s < —(f), and that it 
starts to increase exponentially with the bias voltage 
around ev s ~ —(f). At the bias voltages of ev s > 0, 
0o,o->n increases again with the powers of the bias volt- 
age. The enhancement for tro.o^n caused by the expo- 
nential increase in the range of — 4> < ev s < is up to 
(Al/As) 40 /^. The enhancement factor (Al/As) 407 ' 1 " 
is expected to reaches 2 15 « 3.3 x 10 4 from the estimated 
values. Though the estimations for Al/As and </> do not 
have a high accuracy, we know that the rate of the local 
transformation increases steeply around eV s — e a + A n , 
and the enhancement reaches to the order of 10 4 . A„ 
for the flip-flop motion and e a are estimated to be about 
O.leV and 0.6eV, respectively, from the height of the en- 
ergy barrier for the flip-flop motion and the width of the 
7r*-band and the STS results H, . Using the estimated 
values of A„ and e a , we get V s « 0.7V at which the con- 
dition ev s — is satisfied. This value of 0.7V for the 
bias voltage shows good agreement with the experimen- 
tal results H,@ where the transformation to the p(2x2) 
structure for V s < 0.7V is so slow as to be hardly ob- 
served, and that for 0.8V < V s is fast enough to be well 
observed. 

The excitation rate of the vibrartion of the dimer at the 
position r away from the dimer beneath the tip, c^o^n 
is obtained in the lowest order in \8e\ in the essentially 
same way as cro,o— >n : 

<?y,o— m 

= 2 T ^2 (1 " /(£q ' " £F )) / ( £ f - eVs - £p ) 

x I (a', n, r\H cv (l + G c - v (e')H cv ) \p', 0, 0) | 2 
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xS(e a i + nhco — s p >) 
= n\ 1 J / deA t (e + n?iw)A s (e) 

n-l 

x|G r0 (£ + nM! 2 II |Grr(£+iM| 2 , (6) 
3=0 

where G rr (e) = (a r \(e + iO + — iJ c ) _1 |a r ) = Goo(e) and 
G r o(e) — (a r \(e + ?0 + — i/ c ) _1 |ao}- In order to obtain 
G r o(e), we rewrite the Hamiltonian H e by the bases of 
| do), | Or) and their orthogonalized states \q): 

H c = ^2 £ q c \°q + £aala Q + e a ala r 
q 

+ ^o q a\cq + H.c.) + r rq aic q + H.c.) 

q q 
+ (W r0 ala + H.c), (7) 

where the terms concerning with the tip are omitted. 
The transfer integrals of T rq and T$ q in eq. (7) have the 
magnitudes of same order as r Q in eq. (1) for large r. 
Using the rewritten form of H c in eq. (7), we obtain 
easily G r o(e) as G r0 (e) = (e - e a + i0 + - Goo) _1 (CrO + 
W r o)Goo, where G 00 = J2 q l r o g | 2 (e ~ £ q + i0 + ) _1 and 
G r o = J2 q ^rq^Q q (£ — £ q + i0 + )~ 1 . The integration in 
eq. (6) is approximately evaluated at eV s — rihuj — X n 
same way as for cro,o^n- The excitation rate oy^o— >n is 
approximately obtained as 

C r ,0-»n ~ |(— iA r0 + Wro)| 2 Ao <T ,0-m) (8) 

where A 00 = tt J2 q \T 0q \ 2 S(ev s + s a - e q ) and A r0 = 
7r r rg rQ g <5(ew s + e a — £q)- The absolute value of W r o 
almost vanishes out within the distance of a few dimers. 
For the one-dimensional band, both the absolute values 
of A r o and W r o vanish on the dimer rows other than the 
row including the dimer beneath the tip. On the dimer 
row including the dimer beneath the tip, the absolute 
value of A r0 does not decay with the absolute value of 
r for the one-dimensional band, because only a few dis- 
crete points of q contribute to the summation of q for A r o 
for the one-dimensional band. We know that because of 
the quasi-one-dimensional charactoer of the 7r*-band, the 
transformation to the p(2x2) structure for V s > is re- 
stricted almost into the dimer row including the dimer 
beneath the tip and the rate of the transformation de- 
cays much slowly with the distance from the dimer be- 
neath the tip. The transformed region into the p(2x2) 
structure along the dimer row is expected to extend to 
be of the order of the mean free path of the electron in 
the 7r*-band. 

For the two-dimensional band, the excitation rate 
a r fi— , n is also obtained in the same form as in eq. (8). 
In spite of the fast decay of the excitation rate on r, 
W r o does not vanish on both the direction of parallel and 
perpendicular to the dimer row in the two demensional 



band. The absolute value of W r o is effective only within 
the distance of a few dimers. For the two-demensional 
band, A r0 shows the quasi-isotropic power-law decay on 
r, because in the two-dimensional band, the continuous 
points of q on a line in the g-space contribute to the sum- 
mation of q for A r o, and the summation of q is replaced 
by integral on the line in the g-space. These results for 
the position dependence of the rate for the two demen- 
sional band reproduce well the experiment results 0,0 of 
the local transformation to the c(4x2) for V s < 0, where 
the transformed region is restricted within the range of 
a few nm away from the dimer beneath the tip and the 
transformed region has the quasi-isotropic shape. 

In conclusion, we have investigated theoretically the 
coherent excitation rate of the dimer vibration, using 
the Hamiltonian consisting of the terms of the electron 
system and the electron-vibration coupling terms. We 
show that the transformations of the local structures by 
STM is driven by the dimer vibration excited coher- 
ently by the STM current. The sample bias voltage above 
which the transformation to the p(2x2) structure is able 
to be observed, is semiquantitatively reproduced by the 
quasi-one-dimensional character of the 7r*-band. For the 
one-dimensional band, the excitation rate is shown to 
have a term not decaying on the distance from the STM 
tip. The contrastive diffrence of the shapes of the trans- 
formed region depending on the sign of the bias voltage 
observed by the experiments 0, Q is explained by the 
dimensionality of the 7r-band and the 7r*-band. 



[1] S. Ferrer, X. Torrelles, V. H. Etgens, H. A. van der Vegt 
and P. Fajardo, Phys. Rev. Lett. 75, 1771 (1995) . 

[2] C. A. Lucas, C. S. Dower, D. F. McMorrow, G. C. L. 
Wong, F. J. Lamelas and P. H. Fuoss, Phys. Rev. B 47, 
10375 (1993). 

[3] T. Sato, M. Iwatsuki and H. Tochihara, J. Electron Mi- 

crosc. 48, 1 (1999). 
[4] Y. Yoshimoto, Y. Nakamura, H. Kawai, M. Tsukada and 

M. Nakayama, Phys. Rev. B 61, 1965 (2000). 
[5] Y. Yoshimoto, Y. Nakamura, H. Kawai, M. Tsukada and 

M. Nakayama, Surf. Rev. Lett. 6, 1045 (1999). 
[6] Y. Nakamura, H. Kawai, Y. Yoshimoto and M. Tsukada, 

Surf. Sci. 493, 361 (2001). 
[7] H. Kawai, Y. Yoshimoto, H. Shima, Y. Nakamura and 

M. Tsukada, J. Phys. Soc. Jpn. 71, 2192 (2002). 
[8] Y. Takagi, Y. Yoshimoto, K. Nakatsuji and F. Komori, 

J. Phys. Soc. Jpn. 72, 2425 (2003). 
[9] Y. Takagi, Y. Yoshimoto, K. Nakatsuji and F. Komori, 

Surf. Sci. 559, 1 (2004). 
[10] H. Kawai and O. Narikiyo, J. Phys. Soc. Jpn. 73, 417 

(2004) (Errata: J. Phys. Soc. Jpn. 73 No. 6 (2004)). 
[11] B. N. J. Persson and M. Persson, Solid State Comm. 36, 

175 (1980). 

[12] B. N. J. Persson and A. Baratoff, Phys. Rev. Lett. 59, 
339 (1987). 

[13] S. Gao, M. Persson, and B. I. Lundqvist, Phys. Rev. B 



■5 



55, 4825 (1997). 
[14] N. Mingo and K. Makoshi, Surf. Sci. 438, 261 (1999). 
[15] N. Mingo and K. Makoshi, Phys. Rev. Lett. 84 3694 

(2000). 

[16] K. Makoshi, N. Mingo, T. Mii, H. Ueba, and S. 

Tikhodeev, Surf. Sci. 493, 71 (2001). 
[17] N. Mingo, K. Makoshi, T. Mii, and H. Ueba, Surf. Sci. 



482-485, 96 (2001). 
[18] B. N. J. Persson, H. Ueba, Surf. Sci. 502-503, 18 (2002). 
[19] G. P. Salam, M. Persson and R. E. Palmer, Phys. Rev. 

B 49, 10655 (1994). 
[20] M. Rohlfing, P. Kriiger and J. Pollmann, Phys. Rev. B 

54, 13759 (1996). 



